Determination of fragmentation functions and 
their uncertainties from e + + e ~ ^ h + X data 

M. Hirai*, S. Kumano 1 ^*, T.-H. Nagai** and K. Sudoh 1 " 

* Department of Physics, Tokyo Institute of Technology 
\ Ookayama, Meguro-ku, Tokyo, 152-8550, Japan 

' Institute of Particle and Nuclear Studies, High Energy Accelerator Research Organization (KEK) 
1-1, Ooho, Tsukuba, Ibaraki, 305-0801, Japan 
** Department of Particle and Nuclear Studies, Graduate University for Advanced Studies 
O ' 1-1, Ooho, Tsukuba, Ibaraki, 305-0801, Japan 



(N 
O 

Q 



Abstract. Fragmentation functions are determined for pions, kaons, and nucleons by a global anal- 
ysis of charged-hadron production data in electron-positron annihilation. The optimum functions 
are obtained in both leading order (LO) and next-to-leading order (NLO) of a s . It is important that 
uncertainties of the fragmentation functions are estimated in this work by the Hessian method. We 
found that the uncertainties are large at small Q 2 and that they are generally reduced in the NLO 
Sij in comparison with the LO ones. We supply a code for calculating the fragmentation functions and 

, their uncertainties for the pions, kaons, and nucleons at given z and Q 2 . 
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INTRODUCTION 



Fragmentation functions are used for calculating cross sections of high-energy hadron- 
production processes. They describe how a hadron is produced in the final state from a 
parent quark, antiquark, or gluon. Recently, they are becoming increasingly important 
because semi-inclusive hadron-production processes are investigated for finding the ori- 
gin of nucleon spin and properties of quark-hadron matters in lepton-nucleon scattering 
and hadron-hadron collisions. 

There have been studies of determining the functions from hadron-production data in 
electron-positron annihilation [1]. It is known that determined functions are very differ- 
ent between the parametrization groups, for example KKP (Kniehl, Kramer, and Potter) 
and Kretzer. It is unfortunate that uncertainties of the functions were not estimated al- 
though they have been studied in parton distribution functions (PDFs) in the nucleon 
[2, 3] and nuclei [4]. The major purpose of our work is to investigate the uncertainties 
of the fragmentation functions [5] in both leading order (LO) and next-to-leading order 
(NLO) of the running coupling constant a s . 



ANALYSIS METHOD 



The fragmentation functions are determined by analyzing the data for charged-hadron 
production in electron-positron annihilation. The fragmentation function for hadron h is 



defined by the cross section for e + + e — > h +X: 

o t ot dz 

where Q 2 = s with the center-of-mass energy y/s, and o tot is the total hadronic cross 
section. The variable z is given by the hadron energy Eh and the beam energy y/s/2 by 
z = £/,/ {y/s/2) = 2Eh/ Q. The fragmentation process occurs from primary quarks, anti- 
quarks, and gluons, so that the function is expressed by their contributions: F h {z, Q 2 ) = 
Y,iQ{z,OC s ) ®D^{z,Q 2 ), where D^{z,Q 2 ) is a fragmentation function of the hadron h 
from a parton i (= w, d, s, ■ ■ ■, g), Q{z, <X S ) is a coefficient function, and <g) indicates a 
convolution integral [5]. 

In our analysis, the functions are defined at fixed Q 2 (= Qfy as [5]: 

D l ?{z,Q 2 o)=Nh a <{l-z)f } i, (2) 

where Nf, af, and /3/* are parameters. Instead of Nf, we use a second moment M\ as 
a parameter in our analysis because its physical meaning is clearer. They are related 
by Nf = M ] //B{a } / + 2,/3f + 1), where B{a\ + 2,j3f + 1) is the beta function. The 
parameters are determined from a x 2 analysis of the data. The data are taken at large 
Q 2 , typically Q 2 = Aff , and the functions are evolved to experimental Q 2 points by 
the timelike DGLAP equations. The parameters are determined by minimizing the total 
X 2 , X 2 = Y,i{Ff ata ~ F' heo ) 2 /{of ata ) 2 , where Ff ata and Ff eo are experimental and 
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theoretical fragmentation functions, and c4 ata is an experimental error. The uncertainties 
of the obtained functions are calculated by the Hessian method [2, 3, 4, 5]: 

where Htj is the Hessian, is a parameter, and the parameter set at the minimum % 2 is 
denoted by q . 



RESULTS 

The parameters are determined by fitting the data for e + e~ — > hX. As an example, 
the obtained fragmentation functions in the NLO are compared with experimental 
data for the charged pions (k + + n~) in Fig. 1, where fractional differences {Data- 
Theory)/ Theory are shown for the fragmentation functions at the experimental Q 2 
points. There are fourteen parameters in the pion analysis. The number of the data is 
264. The bands indicate the one-c uncertainty ranges estimated by the Hessian method. 
The data are generally well explained by our parametrization; however, the DELPHI 
data significantly deviate from our fit. The total x 2 is 433.5 in the NLO analysis for the 
pion so that # 2 /d.o.f. is 1.73. In the LO analysis, the total % 2 is slightly larger. 

Determined fragmentation functions for K + are shown at Q 2 =l GeV 2 , m 2 , m 2 ,, and 
M| in Fig. 2. The functions in the LO and NLO are indicated by the dashed and solid 



curves, and the uncertainties in the LO and NLO are shown by the light- and dark- 
shaded bands, respectively. It is important that the uncertainties are now shown in this 
work for the fragmentation functions. At small Q 2 (~ 1 GeV 2 ), the uncertainties are 
large, especially in the LO. The uncertainty bands become smaller in the NLO. It is very 
difficult to determine the gluon function; however, an improvement can be seen in the 
NLO because the uncertainty becomes much smaller. The functions are relatively well 
determined at large Q 2 (~ A/|), where the e + e~ data are taken. 
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FIGURE 1. Comparison with experimental data for pions [5]. 





FIGURE 2. Determined fragmentation functions for n + at Q 2 = 1 GeV 2 , m 2 ., m\, and M| [5]. 
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Similar analyses have been done for the 
kaons and proton/anti-proton. In both cases, 
the values of # 2 /d.o.f. are about two, and 
the data are generally well explained by our 
parametrization. As an example, the functions 
for K + are shown in Fig. 3 with their uncer- 
tainty bands. The notations are the same as 
the ones in Fig. 2. The functions have large 
uncertainty bands also for the kaon at small 
Q 2 . The large uncertainties suggest that the 
hadron-production cross sections have large 
ambiguities if they are calculated in the small 
Pt region. The obtained functions for the pro- 
ton as well as the detailed discussions on the 
pions and kaons are found in Ref. [5]. 

From the x 2 analysis of the data for 
e + e~-^ hX, the optimum fragmentation func- 
tions are determined for the pions, kaons, and FIGURE 3. Fragmentation functions for K+ [5]. 
nucleons. It is especially important that the 

uncertainties of the obtained functions are calculated in this work. The results indicate 
large uncertainties at small Q 1 , which suggests large ambiguities in extracting infor- 
mation, for example on nucleon spin and quark-hadron matters from hadron-production 
cross sections in lepton-hadron and hadron-hadron reactions. In such studies, it is impor- 
tant to indicate reliable regions of the fragmentation functions by using the uncertainties. 
We supply our code for calculating not only the optimum fragmentation functions but 
also their uncertainties at given z and Q 2 for general users [5]. 
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